Astronomy & Astrophysics manuscript no. 100340ms 


©ESO 2010 


July 22, 2010 





o 



Photometry and spectroscopy of GRB 060526: A detailed study of 
the afterglow and host galaxy of a z=3.2 gamma-ray burst ^ 

C. C. Thone''2, D. A. Kann^ G. Johannesson'^ ^ J. H. Selj^, A. O. Jaunsen^ J. P. U. Fynbo\ C. W. Akerlof^ K. S. 
Baliyan^ C. Bartolini'', I. F. Bikmaev^", J. S. Bloom", R. A. Burenin^^^ g g Cobb"'^^ S. Covino^, P. A. Curran^'^ i^ 

H. Dahle^ A. Ferrero^'^, S. Foleyi'^-^^ J. French^^ A. S. Fruchterl^ S. Ganesh^ J. F. Grahaml^ G. Greco'', A. 
Guarnieri'', L. Hanlon''', J. Hjorth^ M. Ibrahimov'^, G. L. Israel^", P Jakobsson^^ , M. Jelmek^^, B. L. Jensen^ U. G. 
J0rgensen23'24, i. m. Khamitov^^, T. S. Koch^*', A. J. Levan^^, D. Malesani^ N. Masetti^^ S. Meehan^^, G. Melady^^, 

D. Nanni^", J. Naranen^^, E. Pakstiene^°, M. N. Pavlinsky^^, D. A. Perley", A. Piccioni^, G. Pizzichini^^ , A. 
Pozanenko^^, P W. A. Roming^^.^i, w. Rujopakarn^^^ y. Rumyantsev3^ E. S. Rykoff^'*■^^ D. Sharapov^^, D. Starr^^, R. 
A. Sunyaevi2-36, H. Swan^ N. R. Tanvir^^ R Terra^^ P M. Vreeswijk\ A. C. Wilson^'', S. A. Yost"*", and F Yuan'^ 

(Affiliations can be found after the references) 
Received 8 June 2008 / accepted 17 July 2010 



O 

(N 



I 

o 



> 

00 



o 

00 

o 



ABSTRACT 

Aims. With this paper we want to investigate the highly variable afterglow light curve and environment of gamma-ray burst (GRB) 060526 at 
z = 3.221. 

Methods. We present one of the largest photometric datasets ever obtained for a GRB afterglow, consisting of multi-color photometric data from 
the ultraviolet to the near infrared. The data set contains 412 data points in total to which we add additional data from the literature. Furthermore, 
we present low-resolution high signal-to-noise spectra of the afterglow. The afterglow light curve is modeled with both an analytical model using 
broken power law fits and with a broad-band numerical model which includes energy injections. The absorption lines detected in the spectra are 
used to derive column densities using a multi-ion single-component curve-of-growth analysis from which we derive the metallicity of the host of 
GRB 060526. 

Results. The temporal behaviour of the afterglow follows a double broken power law with breaks at / = 0.090 ± 0.005 and / = 2.401 ± 0.061 days. 
It shows deviations from the smooth set of power laws that can be modeled by additional energy injections from the central engine, although some 
significant microvariability remains. The broadband spectral-energy distribution of the afterglow shows no significant extinction along the line of 
sight. The metallicity derived from S 11 and Fe II of [S/H] = -0.57 ±0.25 and [Fe/H] = -1.09±0.24 is relatively high for a galaxy at that redshift but 
comparable to the metallicity of other GRB hosts at similar redshifts. At the position of the afterglow, no host is detected to F775W(AB) = 28.5 
mag with the HST, implying an absolute magnitude of the host M(1500 A)>-18.3 mag which is fainter than most long-duration hosts, although 
the GRB may be associated with a faint galaxy at a distance of 1 1 kpc. 

Key words, gamma rays: bursts - gamma rays bursts: individual: GRB 060526 - galaxies: high redshift - interstellar medium: abundances 



1. Introduction 

I Gamma-ray bursts (GRBs) and their afterglows offer a powerful 
. tool to probe the high-redshift universe, both through photom- 
etry and spec troscopy. The standard fireball model of GRB af- 
] terglo ws (see 'Zhang'' 2007t iMeszaros etal]|2006t iGehrels et al] 



I2009L for recent reviews on the topic) predicts a smooth tem- 
' poral evolution, and the resulting afterglow light curve can be 
empirically described by a joint smoothly-broken power law, 
the so-called Beuermann equation (Beuermann et al. 1999). For 
many GRBs in the pie-Swift era, when temporally dense af- 
terglow photometry was obtained, the afterg low evolution was 
foun d to be smooth as, e.g., in GRB 0208 1 3 dLaursen & Staned 
[2003), GRB 030226 jKlose et all I200I and GRB 041006 
(Staneketal. 2005) (for the complete pie-Swift sample see 
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tained with the NASA/ESA Hubble Space Telescope under proposal 
11734 (PI: Levan). 



IZehetal.ir2006h . Out of the total sample of 59 afterglows of 
that paper, only four of the GRBs analysed showed significant 
deviations from the expected decay (though we note that about 
half of the 59 afterglows were not sampled well enough to al- 
low any conclusion s). One might be explained by microlens- 
ing (GRB 000301C. lGarnavich et an i2000) and another one by 
assuming an inhomogeneously emitting surface (GRB 011211, 
Jakobsson et al. 2004). The other two GRBs which are the only 
ones showing long-lasting s trong deviations are GRB 021004 
(e.g. Ide Ugarte Postigo et al.h2005i, an d references therein) and 
GRB 030329 (e.g. iLiEkiiLetdJ |2004|, and references therein), 
which incidentally also have the densest optical monitoring. 

With the launch of the Swift satellite and its rapid localiza- 
tion capabilities (Gehrels et al. 2004), the number of highly vari- 
able light curves has increased considerably, though there are 
still examples of very smooth light curves (e.g. that of GRB 
080210, which showed a smooth behaviour with a high sam- 
pling rate of 1 s, A. De Cia, priv. comm.). Some light curves 
show small, achromatic bum ps overlying the smo oth power law 
decay (e.g., GR B 050502 A. Guidorzi et al.l l2 00l GRB 061007. 



iMundell et al.]l2007c GRBs 090323 and 090328, 



McBreen et al.l 
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I2OIOI) . early bumps with chromatic ev olution (e.g., GRB 
061126. iPerlev et al.l l2008a': GRB 071003. lPerlev et al.ll2008bh . 
"steps" due to energy injection episo des (e.g., GRB 07 0125, 
lUpdike e t al.' 2008: GRB 071010A, Cov ino et alJ [ 2008; GRB 
080913, ' Greiner et__al 
ICenko et al 



201(f 

time rebrightenin gs of up to 
050721, AntonellieTai] |2006|; 



.,2009b. GRB 09092 6A. iRau et aL,201Q: 
Swenson et al.1 1201 Ot) or powerful late- 
several magnitudes (e.g ., GRB 
GRB 060206^ Wozniak et all 



2006; Monfardini et al. 2006; Stanek et al. illOOTt GRB 07 0311, 
Guidorzi et al.. .2007: GRB 071003. .Perlev et al.1 120081 . The 
early time domain, which can now be routinely accessed by 
rapid follow-up in the Swift era, has yielded mor e types of vari- 
ability , like rising afterg lows (e.g., GR B 060418. rMolinari et al .l 

12009; GRB 060607A, 
-rrr, ' „ „„„ 



2007t GRB 060 605. 



iFerrero et al 



Nvsewander eTap 12009*. 'Molinari et al.' 2007; G RB 081008 



Yuan et al.l l2010t seelOates et al. 2009 and Rvko ff et al.1 l2009l 



for further examples) and short-term variability directly linked 
to the p rompt emission (e.g., GRB 041219A, I Vestrand et al .l 



[2005; B lake etal.1 120051; GRB 050 820A. Vestrand et al 



GRB 80319B, lRacusin et"al]|2008t GRB 080129. .Greiner etalj 
l2009ah . In all these cases, dense photometric follow-up during 
the periods of variability was needed to characterize the phe- 
nomena involved. 



In addition their use in to studying the GRB phenomenon 
itself, GRB afterglows can be used to study their galactic en- 
vironment through absorption line spectroscopy of material in 
the line-of-sight towards the GRB. Above a redshift of around 
z ~ 0.2, resonant absorption lines from elements present in the 
interstellar medium (ISM) such as Mg, Zn, Fe, Si, C, and S are 
shifted into the optical regime and can be studied with ground- 
based telescopes. For 10 bursts between a redshift of z = 2 to 6, 
the metallicity along the line-of-sight in the galaxy could be ob- 
tained so far (e.g. Savaalio 2006; Fvnbo et al. 2006; Pric e et all 
I2OO7L iProchaska et al.ll2007i: iLedoux et al.ll2009 . and references 
therein). The values are usually below solar, but higher than 
for QSO absorbers at comparable redshifts, some of them even 
higher than theoretical limits for the formation of collapsars 
(Woosley & Heger 2006). The difference to QSO absorbers can 
be explained if GRB sightlines probe denser parts of the galaxy, 
or if GRBs reside in g alaxies w ith high er masses and there - 
fore higher metallicities jFvnbo e t al. 2008; lPontzen et al.l2010h . 
Both for QSO and GRB absorbers, there see ms to be a metal- 
licity evolution w ith redshift (Savaglio 2006; F vnbo et ani2006t 
iPrice et al.l2007h . although the slope is different for the two sam- 
ples. GRB hosts s eem to show a low extinction along the line- 
of-sight (Prochask a et al.ll2007h . however, relative abundances 
of heavier elements indicate that some of the ions must be de- 
pleted onto dust grains and the depletion pattern resembles the 
one found in th e warm disc and halo of the Milky Way (MW; 
ISavaglioll2006l) . How this can be explained together with the 
low extinction as also derived from the spectral energy distri- 
bution (SEP) of t he afterglow ( Kann et al.l[2006t IStarling et aP 
l2007tlKannetd]l20rQ) . is still an open question. One possi- 
ble solution to this problem is destruction of the dust present 
in the line-of-sight by the GRB and afterglow radiation (e.g., 
[Waxman & Draine 2000; Perna & Lazzati 2002), though no 
strong evidence has ever been found for this. Furthermore, there 
is clearly a strong observational bias involved, as those GRB af- 
terglows with successful spectroscopy, especially in the case of 
high-resolution observations, are those which have only low ex- 
tinctions and thus relatively bright afterglows (Fvnbo et al. 2009; 
iKann et al]|2010l) . although rapid observations with large tele- 
scopes can achieve detailed spectroscopy of highly extinguished 



afterglows, as in the case of GRB 0806 07 (IProchaska et al.l2009l; 
ISheffer et al.ll2009HFvnbo et al.ll2009h . 

GRB 060526 was detected by the Swift satelUte on May 
26.686458 (16:28:29.95 UT). The satelHte slewed immediately 
to the burst, d etecting both the X-ray and the optical afterglow 
(ICampana et a l. 2006a). The BAT instrument on Swift measured 
two emission episodes. The first one lasted 13.8 s and consisted 
of two FRED (fast rise exponential decay) peak s, followed by a 
second symmetric peak between 230 and 270 s dCampana et al.l 
l2006bh . The second peak was coincident with a giant X-ray 
flare followed by a softer flare at 310 s (Campana et al. 2006c|) 
also d etected in the v b and by the UVOT telescope on-board 
Swift (lBrownetalJl2006 ). The gamma-ray fluence was (4.9 + 
0.6) X 10"^ erg cm"^ during the first emission episode and 
(5.9 + 0.6) X 10"' erg cm"^ during the second, the peak flux 
of the second episode was however only half of the peak flux of 
the first one. The photon index of the two epochs changed from 
1 .66 + 0.20 to 2.07 + 0.18, thus show ing the typical hard-to- 
soft evolution (iMarkwardt et al.l[2006h . The Watcher telescope 
provided the first ground based detection with R ^ 15 mag 
(French & Jelmek 2006) 36.2 s after the trigger ROTSE obser- 
vations sho wed a plateau for se veral thousand seconds after the 
GRB onset (iRykoff et al."2006'). A redshift of z = 3.21 was de- 
termined by Berger & Gladders (2006) with the Magellan/Clay 
telescope. The brightness of the optical afterglow allowed for a 
dense monitoring which revealed a complex light curve structure 
including several flares (Halpern et al. 2006j|a) and a steepening 
attributed to a jet break (Thone et al.ii2006i) . 

In this paper, we approach the analysis of GRB 060526 
from two directions: Through modelling of the very detailed 
optical light curve and late optical imaging of the field to de- 
tect the host in Sec. 3, and analysis of low to medium reso- 
lution spectroscopic observations of absorption lines along the 
line of sight (Sec. 4). Throughout the paper, we follow the con- 
vention Fy(t) oc f~"v ~^, and use WMAP concordant cosmology 
(ISpergel et al.ll2003h with Hq = 71km s ' Mpc ', Qm = 0.27, 
and Qa - 0.73. Uncertainties are given at 68% confidence level 
for one parameter of interest unless stated otherwise. 



2. Observations 

2.1. Photometry 

In order to get a good coverage of the light curve, we obtained 
data using several different telescopes around the world. Our 
complete data set comprises a total of 412 points from the UV 
to K-band, one of the largest photometric samples of an opti- 
cal/NIR afterglow in the Swift era. 

The earliest dataset was obtained by the Watcher tele- 
scope, located at Boyden Observatory, South Africa, start- 
ing 36.2 s after the burst, followed by Swift UVOT start- 
ing 86 s after the trigger. Early ground based optical data 
were obtained with the ROTSE-IIIc 0.3m telescope at the 
H.E.S.S. site at Mt. Gamsberg, Namibia, the 1.5m telescope 
on Mt. Maidanak/Uzbekistan, the 2.6m Shajn telescope at 
CrAO (Crimean Astrophysical Observatory/Ukraine), the TNG 
(Telescopio Nazionale Galileo) on La Palma equipped with 
DOLoRes, the 1.2m MIRO telescope on Mt. Abu/India, with 
BFOSC (Bologna Faint Object Spectrograph & Camera) at the 
G. D. Cassini 152 cm telescope of the Bologna University un- 
der poor conditions and with the RTT150 (1.5m Russian-Turkish 
telescope , Bakirlitepe, Turkey), the RTT150 data are also pre- 
sented in lKhamitov et alJ (l2007h . 
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The light curve was followed up sparsely every night over 
nearly a week with DFOSC (Danish Faint Object Spectrograph 
and Camera) on the Danish 1.54m telescope on La Silla/Chile 
under partially photometric conditions and with MOSCA 
(MOSaic CAmera) and ALFOSC (Andalucia Faint Object 
Spectrograph and Camera) at the Nordic Optical Telescope on 
La Palma. Furthermore, two epochs were obtained with the 
Tautenburg L34m Schmidt telescope and two sets of images 
were taken several days after the GRB with Keck/LRIS simul- 
taneously in the Kron-Cousins R and the Sloan g' bands, the 
second observation was performed at high airmass under bad 
seeing conditions. Late images were obtained with FORS2 at 
the VLT on Paranal/Chile on Feb. 23, 2007 and Mar. 30, 2008 
in the Rc band with exposure times of 2500 and 7500 s, respec- 
tively, to look for the host galaxy. We also took observations 
with the Hubble Space Telescope on 9 August 2009, utilizing 
the Advanced Camera for Surveys with the F775W (roughly 
SDSS /') filter (see Fig. |6]l. A total of 7844 s of observations 
were obtained in six dithered exposures. These were reduced via 
multidrizzle in the standard fashion. 

Near infrared data were collected with ANDICAM and 
the 1.3 m SMARTS telescope (Small and Moderate Aperture 
Research Telescope System) at CTIO under non-photometric 
conditions as well as with the robotic L3m PAIRITEL telescope 
on Mt. Hopkins. 

The UVOT data were reduced and analysed using the stan- 
dard UVOT tasks within the heasoft package. For the photo- 
metric calibration of the ground-based data, we determined the 
calibrated magnitude of six comparison stars in the field us- 
ing photometric zero points from DFOSC in the V, Rc and Ic 
bands (see Table |3]l. These stars were then used to perform rel- 
ative Point-Spread Function (PSF) photometry to get the cali- 
brated magnitude of the afterglow. For some of the late NOT 
images as well as the faint MIRO detections, though, we applied 
relative aperture photometry using a circle of 20 pixels diame- 
ter (and an annulus of 10 pixels for the sky). For the B-band, 
where no DFOSC data were available, we took zero points for 
only three comparison stars from t he SDSS (se e Table O, con- 
verting them with the equations of iJester et al.l (l2005l) . g' zero 
points were taken from the SDSS, and they are given as AB 
magnitudes. For the RTT 150 data, we used one USNO-Bl star 
as reference that was calibrated using Landolt standard stars. 
The results are in full agreement with the rest of the data set. 
The Watcher data were an alysed using a dedicated photome- 
try pipeline (Ferr ero et al.l 12010). The J, H and Ks band data 
from SMARTS/ANDICAM and PAIRITEL were calibrated us- 
ing three and ten nearby stars, respectively, from the 2MASS 
catalogue. For the H and Ks band data, we used smaller aper- 
tures and applied aperture corrections to reduce the influence of 
the highly variable background. 

All data and upper limits are given in TablelT] the data are not 
corrected for Galactic extinction. Note we give the g' and HST 
FllSW magnitudes in AB magnitude s. For the final h ght curve 
fitting, we add Br'Rci' band data from Dai et al.'('2007). We shift 
the B and Rc band data of Dai et al. (2007) by 0. 1 magnitudes to 
bring it to our zero point. The multi-color light curves are shown 
in Fig.[T] 



2.2. Spectroscopy 

Spectra were obtained with FORSl/VLT on May 27 from 9 to 
12 hours after the burst. Four different grisms cover the wave- 
length range from 3650 to 9200 A. For all four grisms a I'.'O slit 



was used which provides resolutions between 2.4 and 11.1 A. 
Reduction, cosmic ray removal, extraction and wavelength cali- 
bration were performed using standard tasks in IRAlfl The final 
spectra were then normalised as no absolute flux calibration was 
needed. In order to improve the S/N, we combined the datasets 
taken with the same grism weighted with their variance. A sum- 
mary of the spectroscopic observations is given in Table |2] 



3. Prompt emission 

Up to two hours after the burst, the light curve features an opti- 
cal flare contemporaneous to the XRT/BAT flare at ^ 250 s and 
a following plateau phase (iBrown et al ] l2006HFrench (S: Jelmekl 
2006; Rykoff'et al. 200d). We reduced the BAT event data using 
standard procedures within the software provided by HEASOFT 
(version 6.1). The XRT observations were reduced using the 
standard xrtpipeline (version 0.10.4) for XRT data analysis 
software using the most recent calibration files. The spectral 
d ata of BAT and XRT were analysed with XSPEC version 11.3 
(lArnaudil 996). Th e X-ray Galactic col umn density was fixed to 
5.02 X 10^" cm-^ (iKalberlaetal .112005 ). We estimated the late- 
time extragalactic column density by fitting the XRT PC data 
from 517 s to 1.34 x 10^ s post trigger, where spectral evolution 
is negligible, using an absorbed power law mod el. We only find 
an up per limit of Nhx < 9.8x10^' cm"^ (see also lCampana et all 
' 2010h . 

In Fig. |2]we compare the timescale of the flares at high en- 
ergies (BAT and XRT) with contemporaneous optical data. The 
first, stronger flare is seen both by the XRT and BAT whereas 
the second, softer flare is only visible in the XRT data. Our ear- 
liest optical data are also coincident with the two flares. A broad 
early bump is observed in the optical light curve which precedes 
the second episode of BAT emission and the major X-ray pulse. 
The Watcher data only show a plateau during the high-energy 
flares due to a relatively low time resolution. The higher time- 
resolution of the UVOT data shows that there are also two sig- 
nificant flares in the optical. The first optical flare is contempo- 
raneous with the peak of the BAT/XRT flare within 2cr of the 
temporal error. Then, there may also be a small optical bump 
(significance only «; 1.5cr) at the time of the second (XRT-only) 
flare, whereas the second significant optical flare occurs » 90 s 
after the XRT-only flare implying that these two events are prob- 
ably not connected. Power-law fitting of the optical light curve 
(with respect to the BAT trigger time To) shows the slopes are 
very steep. The rising slope of the first flare is a^i - -12.8 + 3.0, 
and the decay slope is a^i = 5.8 ± 0.8, or, if one takes the bump 
as a second flare, the two decay slopes are a^n = 8.2 + 2.0, 
- 8.6 ± 2.5. The rise and decay slopes for the second peak 
are similarly steep, but with larger uncertainties. 

Prompt optical flashes attributed to reverse shocks, a s 
seen in the landmark burst GRB 990123 dAkerlofet al.l 1 19991) . 
have been observed in on ly a few cases si nce the launch of 
Swift, e.g., GRB 06 01 1 IB (iKlotz et alJl2006h and GRB 0601 17 
dJelmek et a n i2006l see lKann etalJl2010l for a recent overview 
of Ught curves with probable reverse shock flashes/steep de- 
cays in the Swift era). Most bursts for which early optical 
data are available display no evidence of reverse shock emis- 
sion. In a number of cases, the optical light curve is domi- 
nated by forward shock emission from very early times, e.g. 
GRB 050401 (Rvkoffetal. 200 1), GRB 060418 an d GRB 
060607A (Molinarietal. 2007; Nysewander et al. 2009), GRB 



060605 (Ferrero et aL2009.) . GRB 061007 (Mundell et al..2007l) 



http://iraf.noao.edu 
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Table 2. Spectra of GRB 060526 from FORS 1/VLT. Times given are the start time of the observations, midexposure times are in 
days after the GRB trigger. 



Time (May Midtime Exptime Grism Spectral range Resolution 



27 [UT]) 


[d] 


[s] 




[A] 


[A] 


01:16:56 


0.37042 


600 


300V 


3650 


-8900 


11.1 


01:29:46 


0.38105 


900 


6001 


7000 


-9200 


4.4 


01:47:28 


0.39335 


900 


1200B 


3860 


-4400 


2.4 


02:05:27 


0.40583 


900 


600V 


4080 


-7200 


4.5 


02:24:00 


0.42043 


1200 


300V 


3650 


-8900 


11.1 


02:57:05 


0.44681 


1800 


1200B 


3860 


-4400 


2.4 


03:30:34 


0.47007 


1800 


600V 


4080 


-7200 


4.5 


04:04:08 


0.49339 


1800 


6001 


7000 


-9200 


4.4 



Table 3. Magnitudes of the comparison stars used for photometry. 



16 - 

17 - 

18 - 

19 1^ 
« 20 k 



^ 21 I- 

(0 

^ 22 - 



23 - 

24 - 

25 - 



# 


Coordinates 






B 


V 


Rc 




Ic 




RA 




Dec 




[mag] 


[mag] 


[mag] 


[mag] 


1 


15:31 


19.6 


+00 


16:59.5 


19.90 ± 0.03 


19.36 + 0.02 


19.02 ± 


0.02 


18.64 + 0.02 


2 


15:31 


22.7 


+00 


17:21.2 


20.30 ± 0.03 


19.51 ±0.02 


19.09 ± 


0.02 


18.71+0.02 


3 


15:31 


20.9 


+00 


16:39.1 


19.97 + 0.03 


19.11 +0.02 


18.55 ± 


0.02 


18.12 + 0.02 


4 


15:31 


14.4 


+00 


17:49.9 




19.78 + 0.02 


19.43 ± 


0.02 


18.81+0.03 


5 


15:31 


16.5 


+00 


17:55.5 




20.34 ± 0.03 


19.96 ± 


0.03 


20.59 + 0.03 


6 


15:31 


18.3 


+00 


18:17.7 




20.49 ± 0.03 


19.45 ± 


0.02 


18.23 + 0.02 




Kg (upper limit) 
H 

H (upper limit) 

• J 

- J (upper limit) 

• Ic 

• i' (Dai) 

• Rc 

" Rj, (upper limit) 

• (Dai) 
CR 

CR (upper limit) 
r' (Dai) 

• white 

• V/v 

• g' 

• B/b 

• B (Dai) 
I I 



b (upper limit) 
u (upper limit) 
uvw1 (upper limit) 
uvm2 (upper limit) 
uvw2 (upper limit) 
' I I 




1E-3 



0.01 0.1 

Days after the burst 



10 



Fig. 1. Light curves of the optical afterglow of GRB 060526. Filled and open circles (detections, Bg'VRc CR IcJHKs) and open 
triangl es (upper limits, uvw2 uvm2 uvwl ubRc CR JHKs ) are presented in this work, and open stars (Br'Ri') are taken from lDai et al.l 
( l2007h . The early light curve especially shows strong variability, with two flares seen in the UVOT v band data. The later light curve 
also shows evident variability, and the large V - R, B - R colours due to the high redshift are evident. Data in this plot are not 
corrected for foreground extinction. 
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0.25 



0.20 - 



0.15 - 



9. 0.10 



0.05 



0.00 



UVOT V/white 
Watcher unfilt. 
Watcher Rc 




14 



200 300 400 
t after the burst [s] 

Fig. 2. The two early flares seen in the BAT and XRT data 
and temporafly coincident optical observations by UVOT and 
Watcher. The second flare is much softer at high energies and 
hence it is only seen in the XRT data. Both flares are reflected in 
the optical data with the first optical flare being coincident with 
the first BAT/XRT flare while the second flare is delayed com- 
pared to the second X-ray flare. The white data have been shifted 
upward by 0.3 mag to the v zero point. 



and GRB 081008 (lYuan et alJ l2010h . An optical component 
of the emission from internal shocks has been invoked to 
explain the correlation between the optical and high-energy 
light curves observed in several bursts, e.g. GRB 041219A 
(Vestrantet al. 2005; Blake et al. 2005; Fan et al. 200 5), GRB 
050820A dVestrand et a l. 2006), and GRB 050904 dWei et al.1 
l2006h . The very rapid variability as measured by the steep slopes 
and the contemporaneous first flare indicates that GRB 060526 
is another case where the early optical light curve is dominated 
by central engine activity. Such time resolution and coincident 
optical and high-energy observations are still rare, however, a 
thorough analysis of the flares and their spectral properties goes 
beyond the scope of this paper. 



4. Afterglow photometry and light curve modelling 

4.1. The multi-color light curve 

To analyse the evolution of the light curve, we use the Rc/CR 
band which has the densest sampling, a total of 319 data points 
in all, from 50 s to more than 7 days after the GRB. For all fits, 
we fix the host magnitude to ot/, = 29 (see § 14.4b . A fit with 
a single power law to these data is very strongly rejected, with 
X^/d.o.f. = 31.89 (with 319 d.o.f.). Even if we remove the ear- 
liest data which are affected by optical flares, starting only at 
400 s after the GRB, the fit is still rejected Or^/d.o.f. = 29.46 
with 308 d.o.f.). A double smoothly-broken power law gives a 
much better fit, but even so, the fit is formally rejected, with 
;^^/d.o.f. - 3.12 with 302 d.o.f. This is due to t he strong vari- 
ability i n the light curve which was first found b y lHalpern et alJ 
(l2006ah and is also discussed in lDai et aP (l2007h . This fit, along 
with the residuals showing the strong variability, is shown in Fig. 
[3] The parameters we find for this fit (ffpiateau - 0.288 ± 0.026, 
ai = 0.971 +0.008, = 2.524+0.052, fM = 0.090+0. 005 days 
tp — 2 .216 + 0.049 days) are concurrent with those of iDai et akl 
(l2007h (who find ai ^ 1.0, 0-2 ~ 2.9 and ti, ^ 2.55 days). 




1E-3 



0.01 0.1 

Days after the burst 



Fig. 3. The Rc-hand light curve of GRB 060526 fitted with a 
double smoothly-broken power law, using data starting at 0.046 
days. The residuals of the fit show an early optical flare and later 
strong variations of up to 0.3 magnitudes. The dashed vertical 
lines mark the break times, the dotted lines the Icr region of 
uncertainty. The second break is soft (n ^ 4) and is thus seen as 
a smooth rollover. 



Due to the high data density, we are able to let the break 
smoothness parameter n vary for the second break (the first break 
had to be fixed to n = 10), and our result (n = 4.3 + 0.7) 
is in agreement with the tentative ai - n correlation found by 
IZeh et al.1 (l2006h . Still, the significant improvement shows that 
the light curve is basically a double smoothly-broken power 
law, and the steep late decay indica t es th at this break is a jet 
break, as first n oted by iThone et"an (l2006t) and also found by 
iDai et al.) (l2007h . The rest frame jet break time of 0.52 days is 
typical for the optical afterglows of the pre-Swift GRB sam- 
ple (Zeh et al. 2006). Jet breaks, a common feature in well - 
monitored pre-Swift optical afterglows (e.g. . IZeh et al.l 120061) . 
have not been found in many Swift afterglow especiall y in the 
X-rays (e.g.,|M angano et al. 20071 iGrupe et aLi r2007: Sat o et al.1 



120071: iRacusin et al.l200 9). and if there are breaks, then the com 
parison between the optical and X-ra y light cu rves show them 
to often be chromatic (iPanaitescu et al . 2006; O ates et al.ll2007h . 
Analysing both the optical and X-ray data of GRB 060526, 
Dai et al. ( |2007|) suggest that the (jet ) break is a chro matic. 

If we take the X-ray data from iDai et al.l ([2007) and fit it 
with the parameters of our double smoothly-broken power law 
fit using X- ray data from 0.06 days onward, we concur that 
the fit is marginally acceptable, with ;i^^/d.o.f. = 1.6 for 40 de- 
grees of freedom. Using our own X-ray reduction, we obtain 



' Note, though, that lDai et al.l ( l2008h argue that optical follow-up in 
the Swift era, especially with UVOT only, often does not persist long 
enough (and to enough depth) to find jet breaks. 
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Fig. 4. Fit to the light curve in X-ray and optical/NIR Vr'Ri'IJHKs bands including a total of five energy injections. Data before 500 
s are excluded from the modelling, as well as the strong X-ray flare. The light curves have been offset by constant factors as given 
in the figure legend for better legibility. The Bg' bands are affected by additional Lyman forest absorption and were not included in 
the fit. For the Bg' bands, the model predictions have been multiplied by 0.50 and 0.55 respectively to match them with the data. 
The residuals are clearly improved in comparison to those plotted in Fig. [3] but short-timescale variations like the one at 0.7 days 
are still not fitted satisfactorily. 



a worse result, with ;^f^/d.o.f. = 1.7 for 22 degrees of free- 
dom. Similarly, with X-r ay data from the Swift XRT repository 
(lEvans et al.l l2007'. '2009'). the fit is rejected (jc^/d.o.f. = 2.3 for 
33 degrees of freedom). The reason is that there are less data 
points, but with smaller error bars, so the outliers are weighted 
more strongly. A fit to the X-ray data alone results in a much 
earlier break time and slopes that are less steep (again in agree- 
ment with iDai et al.l 12007b . but we caution that the late X-ray 
afterglow is only sparsely sampled and shows large scatter. For 
the three different reductions, we derive the following values us- 
ing a smoothly-broken power law (a double-broken power law 
yields no statistical improvement), with n = 10 fixed and no 
host galaxy: ffi = 0.92 ±0.05, a2 = 3.02 ±0.63, f;, = 1.38 ±0.25 
(XRT repository light curve); ffi = 0.90±0.05, = 2.78±0.61, 
tb = 1.34 ± 0.28 (our XRT reducti on); ai = -0.2 9 ± 0.56, 
a2 = 1.69 ± 0.17, th = 0.16 ± 0.05 (Da i et alJl2007l datal. As 
can be seen, the latter fit is very different from the other two 
(which agree fully within error bars). Even using data from 0.01 
days onward (end of the X-ray flare), the fit is still significantly 
different: ai = 0.43 ± 0.06, 02 = 1.69 + 0.17, tb = 0.25 ± 0.06. 



4.2. Modelling the light curve with energy injections 

Motivated by the similarities to such highly variable light curves 
as that of the afterglow of GRB 021004, which was suc- 
cessfully modelled by multiple energy injections ("refreshed 
shocks") (de Ugarte Postigo et al. 2005), we used the code of 
Johannesso n. Bjornsson & (iudmundsson (2006) to model the 
afterglow light curve. The code numerically solves the kinematic 
equations of an expanding shock front and calculates the re- 
sulting synchrotron emission. Relativistic effects are fully taken 
into account and the code supports delayed energy injection 
episodes. Several energy injection episodes are applied as a pos- 
sible scenario to explain the rebrightenings and shallow decay 
of the afterglow. Preliminary results on a smaller data set were 
presented in Johannesson et al. (2009). The number and a time 
range for the energy injections have to be inserted as initial guess 
for the fit and the fit then adjusts them to the best possible time 
within that range and determines the magnitude of the injec- 
tion. As the very early data likely contain some signature of the 
prompt emission, we exclude all data before 400 s as well as the 
very bright X-ray flare. Due to Lyman forest blanketing, data in 
B and g' bands were excluded from the fit. The V band is also af- 
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Table 4. Fits to the spectral energy distribution of GRB 060526. 
Columns are the dust model, the goodness of the fit, the spectral 
slope and the derived extinction in the rest-frame V band. 



Dust 


X'IA.0.L 


P 


Av 


none 


0.425 


0.695 + 0.035 




MW 


0.309 


0.833 ±0.149 


-0.126 ±0.132 


LMC 


0.466 


0.494 ± 0.397 


0.113 ±0.223 


SMC 


0.396 


0.552 ±0.198 


0.055 ± 0.075 



fected by the b lanketing, but was corrected with the model from 
iMadaul {1995') (brightened by 0.17 magnitudes) and included in 
the fit to have a better spectral coverage at early times. 

Figure |4] shows the best fit ix^lA.o.f. - 2.0) found using a 
model with a total of five energy injections: at 0.01, 0.04, 0.07, 
0.35 and 0.60 days. Each energy injection episode adds 2 free 
parametres to the fit in addition to the 6 parametres needed for 
the standard afterglow model discussed below. The initial energy 
injected into the outflow i^Zso = 5;^^ x 10"*^ ergs. The energy 
injections then add 1.8, 4.1, 4.2, 8, and finally 14 times the initial 
energy release Eq to the afterglow, for a total energy release in 



the afterglow of 1.6+°^ x 10^ 



ergs. 



The first three injections are responsible for the shallow af- 
terglow decay between 0.008 and 0.25 days. The quality of the 
data does not allow to discriminate between this three-injection 
scenario and a continuous injection. Using a lower number of 
injections does not yield a satisfactory fit to the observed data. 
Since there are no direct indications of injections in the light 
curve, the time of each of the three injections is not well deter- 
mined. A direct consequence of this is that the energy of each 
individual injection in this phase is not well determined, while 
the total energy released is fairly consistent. The time and energy 
of the last two injections are, however, better constrained by the 
data. 



Further results of the modelling are a high density of the cir- 
cumburst medium of no = 600^5°™ cm"-' (a high value, com- 
parable to the high-z GRB 050904, Frail et al. 2006) , a rather 
low opening angle of - 2?2;|;[;g and an electron index of 
p = 2.1^ "-^ (h ost galaxy extinction is assumed to be negligi- 
ble, see § 14.3b . The peak frequency v„, passes through the op- 
tical/NIR at very early times, while the data are most likely af- 
fected by the prompt emission. The cooling break is between 
the optical and X-rays up to 6 days after the burst, i.e., over the 
whole data span. Furthermore, we find = 2+-^ , x 10"^ and 
sb - 7^g*' X 10"^, with Se being the fraction of the energy in the 
electron population and sb the fraction of the energy in the mag- 
netic field. Note that the definition of has been changed in the 
model from Johannesson, Bjornsson & Gudmundsson (2006) to 
the definition of lPanaitescu & Kumaii (1200 1) to allow for p < 2 
in the model. This model reproduces the global properties of the 
optical/NIR light curves well, with the rather high;t'^/d.o.f. = 2.0 
resulting from rapid variability in the late time afterglow. This is 
incompatible with the homogeneous shock front assumed in the 
numerical code. The X-ray light curve is also fitted reasonably 
well. 
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Fig. 5. Spectral energy distribution of the afterglow of GRB 
060526 in uvw2 uvm2 uvwl uBg'Vr' Rci' IqJHKs, and fits with 
no extinction (straight black line), MW extinction (dotted 
line), LMC extinction (dashed line) and SMC extinction (thick 
dash-dotted line). The fit with MW dust finds (unphysical) 
negative extinction, causing an "emissive component" instead 
of a 2175 A absorption bump. Data beyond 2.2 x 10'^ Hz 
(Vg'Bu uvwl uvml uvw2) were not included in the fit due to 
Lyman forest blanketing, the grey curves represent extrapola- 
tions. The extinction curves of lPeil ([1992) are not coiTectly de- 
fined beyond 3.2 x 10'^ Hz. The UVOT UV filters are upper 
limits only, showing the strong flux decrease beyond the Lyman 
cutoff. The flux density scale is measured at the break time. 



4.3. The spectral energy distribution and liost extinction 

Following the procedures outlined in lKann et alJ (l2006h . we de- 
rive the optical spectral energy distribution (SED) of the GRB 
060526 afterglow and fit it with several dust models (Milky Way, 
MW; Large Magellanic Cloud, LMC; and Small MageUanic 
Cloud, SMC; Pei 1992.) to derive the line-of-sight extinction in 
the host galaxy. Due to the stro n gly va riable light curve, we 
choose the approach iKann et aP ( 2006h used for the SED of 
GRB 030329, and shift the other bands to the Rc-hand zero 
point to derive the colours. With this method, we can also look 
for colour changes. There may be marginal variations in B - R, 
but this colour remains constant within conservative errors, and 
is more sensitive to Lyman forest blanketing, as different filters 
will suffer a different amount of blanketing. 

The SED clearly shows the decreasing flux in the B, g' and 
V bands, and especially in the uvw2 uvm2 uvwl u bands, where 
only upper limits are found, due to the Lyman forest blanketing 
as well as the Lyman cutoff (see Fig. |5]), and we thus do not 
include these filters in our fit. Given the size of the errors, all fits, 
even with no extinction, are acceptable (see Table |4| and we are 
thus unable to prefer one dust model over another. The lack of z 
band data does not allow us to constrain the existence of a 2175 
A bump, therefore we have no evidence in favor of or ruling out 
MW and LMC dust. We do note, though, that MW dust leads to 
slightly negative extinction, while the LMC dust fit yields large 
errors. Thus, we henceforth use the SMC dust fit, as this is the 
most common dust type found in GRB host galaxies, both in 



^ The error limits presented here are Icr estimates, but due to degen- 
eracy in the model parameters, the actual range of parameter values can 
be much larger. 
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vre-Swift iKann et alJl2006t IStarling et al.ll2007h and Swift-em 
(ISchadv et al.l2007HKann et al.l2010l) results. Note that for SMC 
dust, the extinction is within errors as well. 

Assuming the cooling break Vc to lie blueward of the optical 
bands (highly li kely considering the X-ray spectral slope /3x ~ I, 
iDai et al]l2007l and the result of the numerical modelling), the 
standard fireball model gives for the electron power law index 
p - 2/3+ 1 = 2. 10+0.40. This result (albeit with large errors) is in 
agreement with the result from the broadband modelling (§ I4.2| | 
and with the canonic al p - 2.2, and very similar to many other 
GRB afterglows (e.g.|Kametal]2006; Starling et al. 2007). Our 
results are in contrast to Dai et al.l (12007.) . who derive a very 
steep slope /3o - 1 .69^g 4g from Br'i' data only (coiTecting for 
Lyman absorption) and conclude that the optical and the X-ray 
data lie on the same slope. 



4.4. Host search 

In the late-time VLT images, we do not detect any source at the 
position of the afterglow down to a limit of Rc > 27.1 mag 
(which transforms into an absolute magnitude limit of Mr > 
-20.1). There is one source present at ~1'.'5 South-East of the 
afterglow position with Rc - 26.4 + 0.2, which at a redshift of 
Z = 3.221 would transform into a physical offset between GRB 
and host galaxy of ~ 11.5 kpc. A long-slit spectrum covering 
the afterglow position and this galaxy does not show any trace 
at these positions. If this galaxy was associated with the GRB, it 
would be one of the largest offsets between a lon g-duration GRB 
and its host galaxy known dBloom et alj|2002l) . The strong ab- 
sorption lines from the ISM seen in the afterglow spectrum (see 
Sect. |5]) as well as the high circumburst density inferred from 
the numerical modeling (see Sect. 14.2b do in fact not favour a 
large distance from the host galaxy. A position out side of their 
host g alaxies had been inferred for GRB 070125 (iCenko et al.l 
l2008b and GRB 071003 jPerlev et alj l2008bh . however, those 
spectra showed very weak absorption lines, contrary to what we 
observe for GRB 060526. Concluding, we have no dkect spec- 
troscopic evidence for or against an association of the GRB with 
the galaxy. 

We also observed the field with ACS on HST. Astrometry 
was performed relative to our TNG image, yielding a position 
accurate to ~ 0. 1 " . At the location of the burst we do not find any 
evidence for an underlying host galaxy. To estimate the limits we 
place 50 apertures randomly on the sky and measure the standard 
deviation in their count rates. This implies that the underlying 
host galaxy of GRB 060526 is fainter than F775W(AB) > 28.5 
(3cr) or an absolute magnitude of M(1500 A) > -18. 3 mag which 
is fainter than 0.5 L* according to Gab asch et alj (2004). The 
flux density at the host position is 0.004 + 0.005 fiJy. 

Host galaxies of long-duration GR Bs have often been 
found to be faint ir regular galaxies dFruchter et alJ l2006t 
IChristensen et alJ 12004) which are difficult to detect at higher 
redshifts. So far, there are only eight bursts with z > 3 where 
the detection of a host galaxy has been pub lished, namely 
GRB 97121 4 (z ^ 3.4 18, 'Kulkarni et al."l998'), GRB 000131 
(z = 4.500, lAndersen e t al. 2000; Fruchter et al. 2006), GRB 
030323 (z = 3.3718, [Vr eeswiik et al. 2004), G RB 060206 
(z ^ 4.04795, iFvnbo et anr2006 :" Thone et al. 2"008l; IChen at alJ 
20091) . GRB 060210 (z = 3.9133, Fynbo et al. 2009; Perlev et alj 
2009 ). GRB 060605 (z = 3.77 3, Ferrero et al. 2009),' GRB 



GRB 050730, GR B 050908, GRB 060607 A, GRB 07 072 IB (all 
Chen at al."2009), GRB 050904 (iBerger et al.ll2007h and GRB 
060510B (Perley et al. 2009), on the other hand, were not de- 
tected to very deep limits in the optical and the NIR. 

The distribution of pie-S wift Rc band host ga laxy magni- 
tudes peaks at Rc(AB) - 25 dFruchter et al.ll2006l) but extends 
out to 29 with a typical redshift z ~ 1.4. Swift GRBs (and thus 
their hosts) how ever have a higher mean redshift of z - 2.8 
dJakobsson et al.ll20 06a!Fl so the distribution will b e shifted out 
to even fainter magnitudes. lOvaldsen et al.l d2007l) also find a 
higher magnitude for Swift hosts than for pie-Swift bursts by 
comparing the expected detection rate from pie-Swift hosts with 
detections and upper limits derived from imaging the fields of 24 
Swift and HETE II bursts from 2005 - 2006. 



5. Spectroscopy results 

5.1. Line identification 

We detect a range of metal absorption lines as well as a 
Lyman limit system (LLS) originating in the host galaxy of 
GRB 060526. A redsh ift of z = 3.221 was determined in 
iJakobsson et a n d2006bh from a number of these absorption lines 
using the spectra taken with the 600V grism presented in this ar- 
ticle. 

Most of the lines were fitted from the combined 600V spec- 
trum which covers all metal absorption lines detected longwards 
of Ly-a, but provides a higher resolution than the 300V spectra. 
The 6001 spectra only covers the AlII A 1670 line at the same 
resolution as the 600V grism. The range of 1200B is entirely 
within the Lyman-a forest but does not have a high enough reso- 
lution to deblend metal transitions from absorption caused by the 
Ly-a forest lines. We do, however, detect Ly-/3 and Ly-y absorp- 
tion in the 1200B grism. In the blue end of the 300V spectrum, 
one can clearly see the 915 A Lyman break at the redshift of 
the host galaxy. In Fig. |2l we show the combined spectrum of 
the 300V grism with the identified lines indicated. We note that 
we do not detect any intervening system in the sightline towards 
GRB 060526, which is rather unu sual for a GRB sigh tline, in 
particular at that redshift (see, e.g.. lProchter et al]l2006h . Taking 
the strong Mg II A 2796, 2803 doublet, the redshift path probed 
for any intervening system is between z ~ 0.8 and 2.2. The 3cr 
limits on the non-detection of the Mg II doublet vary between 
0.17 A (z = 0.8) and 0.82 A (z = 2.2) (rest frame). 

In order to determine the equivalent width (EW) of the strong 
absorption lines, we fitted the continuum around the lines in re- 
gions that were free of absorption and summed over the absorp- 
tion contained within two times the full-width at half-maximum 
(FWHM) of the lines. For weak lines, we obtained better results 
due to the low S/N by fitting Gaussians. For this fit, we used 
a modified version of the gaussfit procedure provided in IDL0 
which is more rehable in determining the continuum and fitting 
the actual line even if it is slightly blended with a neighbouring 
line. The upper limits on the EWs for a range of ions noted in 
Table |5] was determined from the spectra taken with the 300V 
grism due to the better S/N of those spectra. Between the indi- 
vidual spectra, we do not find any variability in the EW of the 
individual absorption lines. 



090205 (z = 4.6503, D'Ayanzo et al. 2010), and GRB 090323 



(z = 3.568. ICenkoet a l. 2010; McBreenetal. 201 0D. The hos t 
galaxies of GRB 020124 (,Berger et al.»2002t iChen at al.ll2009h . 



* An updated version of the redshift distribution can be found under 
http://raunvis.hi.is/~pja/GRBsample.html - as of July 1, 2010, the mean 
redshift has become lower, z = 2.19, see also Jvnbo et all fcOO^) . who 
find a mean and median of z = 2.2. 

' Available at http ://www.pa.iasf.cnr.it/~rucastro/IDL/Lib/gfit.prol 
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Fig. 6. Left panel: An image from the TNG with the afterglow present. Right panel: The same region observed with HST/ ACS 
and a zoom-in around the region of the afterglow. There is no source detected at the afterglow position itself down to a limit of 
F775W(AB)> 28.5. The brighter spot to the South-East in the zoomed //primage corresponds to the host candidate from the VLT 
mentioned in the text. 




4000 5000 6000 7000 8000 9000 

wavelength [A] 

Fig. 7. The spectrum taken with 300V grism, the two exposures have been combined to improve the S/N. The identified lines as 
well as the main atmospheric absorption bands are labelled. The line at the bottom shows the error spectrum. 



SII ^ 1259, Sill A 1260 and Fell A 1260 are blended and 
cannot be fitted separately. We therefore cannot consider them 
for the derivation of the column density from the curve of growth 
fit as described below and only give the total EW in Table |5] In 
contrast to what is noted in Jakobssonet al. (2006b), we can- 
not reliably detect any fine-structure lines and only give an up- 
per limit for Sill*. Fine-structure lines would be a clear indi- 
cation that the detected absorption lines indeed originate in the 
host galaxy of the GRB, as they are as sumed to be produced 
by UV pumping froni the afterglow (e.g JProchaska et al.ll2006t 
IVreeswiik et alj|2007l: iD'Elia et al.|[2009l) . The redshift derived 



is therefore to be strictly taken as a lower limit only, the de- 
tection of the Lyman a forest redward of the proposed redshift, 
however, excludes a significantly higher redshift for the burst. 



5.2. Column densities from curve of growtin analysis 

Some of the strong absorption lines are saturated, which is a 
problem in low resolution spectra as the damping wings are 
not resolved and Voigt profile (VP) fitting cannot be adopted to 
derive a reliable column density. Furthermore, high resolution 
spectra of GRBs (iProchaskal l2006i) have shown that the strong 
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Table 5. EWs of detected absorption lines and 2cr upper limits on some undetected lines, the EW for the blended systems include 
the contributions from all lines. The column densities were derived from CoG fitted for S II, Si II, Fe II, C IV, Si IV and Al II. Si II* 
denotes a fine-structure line. Upp er and lower limits were d etermined by assuming the ions to lie on the linear part of the CoG. The 
Lya column density is taken from lJakobsson et"an (l2006bh and is based on the 600V grism. 



'lobs 




ID 


z 




log N/cm ^ 


[A] 


[A] 






[A] 




4102.26 


972.54 


Lyy 


3.218 


2.39 ±0.10 




4329.70 


1025.72 


hyp 


3.221 


2.92 ± 0.06 




5131 


1215.67 


hya 


3.221 




20.00 ±0.15 


(5254) 


1242.80 


NY 


(3.221) 


< 0.09 


< 13.66 


(5276) 


1250.58 


SIX 


(3.221) 


< 0.08 


< 14.71 


5291.01 


1253.81 


SII 


3.220 


0.05 ± 0.02 


14.58 ±0.25 


5320.00 


1259.52 


SII 


3.221 


(1.32 + 0.07) 






1260.42 


Sill 


3.221 


blended 






1260.53 


Fell 


3.221 






(5335) 


1264.74 


Si II* 


(3.221) 


< 0.10 


< 12.76 


5497.30 


1302.17 


OI 


3.222 


1.07 ±0.05 


> 15.15 


5506.82 


1304.37 


Sill 


3.222 


0.70 ±0.21 


15.87 ±0.16 


(5559) 


1317.22 


Nill 


(3.221) 


< 0.11 


< 13.59 


5634.24 


1334.53 


CII 


3.222 


1.51 ±0.10 


> 15.55 


5884.87 


1393.76 


Si IV 


3.222 


0.20 ± 0.08 


13.50 ±0.14 


5923.05 


1402.77 


Si IV 


3.222 


0.12 ±0.05 


13.50 ±0.14 


6444.50 


1526.71 


Si II 


3.221 


0.86 ± 0.04 


15.87 ±0.16 


6539.58 


1548.20 


CIV 


3.223 


0.26 ± 0.08 


14.10 ±0.13 


6546.39 


1550.78 


CIV 


3.221 


0.22 ± 0.09 


14.10 ±0.13 


6790.24 


1608.45 


Fell 


3.221 


0.17 ±0.03 


14.28 ± 0.24 


6800.38 


1611.20 


Fell 


3.221 


<0.01 




7052.62 


1670.79 


AlII 


3.221 


1.02 ±0.05 


15.12 ±0.18 


(7863) 


1862.79 


Aim 


(3.221) 


<0.35 


< 12.71 


(8552) 


2026.13 


Znll 


(3.221) 


<0.45 


< 12.73 


(8679) 


2056.25 


CrII 


(3.221) 


< 0.54 


< 13.47 



column density on the linear part of the CoG where the lines are 
not saturated (optical depth tq < 1), but depends on the Doppler 
b parameter of the medium on the flat part of the curve which 
applies for mildly saturated lines. For GRBs, one usually has to 
do a multiple-ion single-component CoG (MISC-CoG) analysis, 
adopting the same effective Doppler parameter for all ions. Here, 
we used all unblended ionic lines that were not heavily saturated, 
namely S II, Si II, Fe II, Al II as well as C IV and Si IV and calcu- 
late the minimum going through the parameter space for the 
column densities of each ion and a range of b parameters. We 
then find the best fit for /? = 39 + 3 km s ' (see Fig. [8]l and the 
column densities as noted in Table |5] Most of the ions lie near 
the linear part of the CoG and are therefore rather independent 
of b. Si II and Al II lie on the flat part of the CoG and the column 
densities are only noted as lower limits as derived from the linear 
part of the CoG. We excluded the saturated O I and C II transi- 
tions from the fit and only list lower limits for these two ions for 
which we take the column density resulting from the linear part 
of the CoG. We also discarded the blended absorption lines of 
Si II A 1260, S II and Fe II even though Si II is the dominant con- 
tribution to the absorption since we detect two other unblended 
Si II transitions at A 1304 and A 1526 A. 

There are several problems connected with the use of the 
CoG that have to be considered. Different ionisation levels 
should actually be treated in separate analyses as they might oc- 
cur in different regions in the absorbing system. However, we 
assume that the absorption takes place in a relatively small re- 
gion of the host galaxy and the resolution of the spectrum does 
resolve different components of the absorbing material in the 
host galaxy. Therefore, we also fit the higher ionisation lev- 
els of Si IV and C IV in the same CoG. Furthermore, they lie 
close to the linear part of the CoG and excluding them from 




log(NfX) 

Fig. 8. Multi-ion single-component curve of growth fit for the 
absorption systems in GRB 060526 using 6 different ions. The 
CoG is expressed in units of the column density "N", the oscilla- 
tor strength of the transition "f" and the wavelength of the tran- 
sition. As comparison, = oo is plotted, the dotted line marks 
the 1 cr deviation from the fit. 



metal absorption lines unresolved in low resolution spectra usu- 
ally consist of a number of narrow, unsaturated components that 
would allow an accurate determination of the column density by 
fitting the different components separately. 

If only low resolution spectra are available, on e has 
to adopt a curve of g rowth (CoG) analysis (Spitzer| Il978t 
ISava^e & Sembachlll996) which directly relates the EW to the 
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the common fit would not change the derived b parameter very 
much. Another problem with doing multiple-ion CoG analysis 
using stron g line s in low-resolution spectra has been noted by 
IProchaskal (120061) . who compared column densities derived from 
low resolution spectra and CoG with high resolution data and 
VP fitting from the same bursts. He found that when including 
saturated lines in the fit, column densities are generally under- 
estimated. Indications for that are if an effective b parameter 
of s> 20 km s ' is found, since strong lines a ctually co n sist of 
a range of components with b < 20 km s^'. ISavagliol (l2006l) 
however performed a similar analysis using CoG and the ap- 
parent optical depth (AOD) method as described in Petti ni et al.l 
(I2002h . which can be applied for medium resolution high S/N 
spectra, and found a good agreement between the two methods. 
The column densities of the saturated lines AlII and Sill are 
infact rather sensitive to the adopted b parameter These lines 
very likely consist of a number of unresolved weaker compo- 
nents which would lie on the linear part of the CoG, the real 
errors should therefore be larger but are difficult to estimate. 

5.3. Metallicity and relative abundances 

Absorption lines that are likely not affected by dust depletion 
can be used to derive a metallicity of the medium in the line of 
sight to the GRB. The least dust-depleted element is Zn, which 
is, however, undetected in our spectra. We then use the relatively 
weak lines S II and Fe II derive relative metallicities compared 
to the hydrogen density with [M/H]=log(NM/N//)-log(NM/N//)0 
using solar abundances from Asplund et al. (2005). Here we 
derive metalHcities of [S/H] = -0.57 + 0.25 and [Fe/H] = 
- 1 .22 + 0.24. 

Fe is usually affected by dust depletion (ISavaee & SembachI 
Il996 f) and corrections have to be adopted. Using the relation be- 
tween the Zn and Fe abundance in Savaglio (2006), we find a 
metallicity of [Fe/H] - -1.09 + 0.24 which marginally agrees 
within the errors with the value derived from sulphur The only 
detected and unblended SII line at A 1253, however, is only 
marginally detected and therefore the EW has large errors. The 
Fe II doublet taken for the CoG fit, in contrast, is also slightly 
blended, but the fit of the stronger component can be consid- 
ered as reliable. Despite the dust depletion, the metallicity de- 
rived from Fe might be the most reliable one in this case and we 
therefore assume a metallicity of [Fe/H] = - 1.09 for the host of 
GRB 060526. 

Independent of the ion used, the metallicity is rather high 
compared to other galaxies at redshift z ^ 3, but am ong the 
typica l metallicities derived for other GRB hosts ( Fvn bo et al.l 
I2006h at that redshift. For those measurements, different ions 
have been used depending on the quality and the wavelength 
coverage of the spectra. Our results show that caution is re- 
quired when comparing the metallicities derived from differ- 
ent elements, as they might be differently affected by dust 
depletion and/or evolution. This is especially true when sat- 
urated lines in low resolution spectra are used to derive the 
column densities and hence the abundance s as it is the case 
for, e.g., GRB 000926 (SavagUo et al."2003") and GRB 011211 
dVreeswijk et al. 2006). For GRB 050401 (Watson et al. 2006) 
and GRB 050505 (iBerger et all 120051) . the authors themselves 
note that due to saturation the reported metallicity is indeed a 
lower limit. 

Our spectra do not allow us to determine a dust depletion 
pattern from the relative abundances of heavier elements, since 
at least four elements out of Zn, Si, Mn, Cr, Fe, Ni or S are nec- 
essary to do such a fit. The difference in the relative abundances 



between S, Si and Fe might suggest the presence of some dust in 
the line-of-sight towards the GRB, however, large extinction is 
excluded from the afterglow SED. On the other hand, this differ- 
ence could also be due to an enhancement in a element produc- 
tior0 which is likely to happen in the young star-forming host of 
a GRB. Generally, GRB host s have higher a/Fe ratios than QSO- 
DLAs (Prochaska etal ] 12007). which can either be interpreted 
as large dust depletion consistent with the higher metallicities of 
GRB sightlines or as a element enhancement. There are indica- 
tions that most of the a/Fe ratio is due to dust depletion traced 
by a large [Zn/Fe] or [Ti/Fe] < 0. The latter is assumed to be a 
clear indicator of dust depletion ([Ti/Fe] < 0) vs. a enhancement 
([Ti/Fe] > 0) (IPessauges- Zavadskv et al.ll2002) . Both elements 
are, however, not detected and the limit on Zn does not lead to 
a strong constraint on the [Zn/Fe] ratio to settle this issue in the 
case of GRB 060526. 

The ratio between high and low ionisation species in the 
spectrum clearly shows that most of the material in the line of 
sight is in a low ionization state. Si IV and C IV have rather low 
column densities and lie close to the linear part of the CoG in 
contrast to their low ionization species Si II and C II which are 
both saturated, whereas Al III is not even detected in our spectra. 
We then derive column density ratios of log(Si IV/Si II) = -2.37, 
log(CIV/CII) < -1.45 and log(Al III/Al II) < -2. 41. From the 
large sample of Swift long GRB afterglow spectra (Fvnb o et al.l 
2009), this seems fairly normal for an average GRB sightline. 
However, the few GRBs occurring in LLSs (log A^Hi/cm"^ < 
20.3) usually show a higher fraction of ionized material com- 
pared to GRB-DLA sightlines. This might either be due to de- 
creased shielding of the highly ionizing afterglow flux by the 
lower hydrogen column density, by a rather special arrangement 
of the GRB inside the host galaxy or simply by a small host 
galaxy. Since our spectra show a low fraction of highly ionized 
material this implies that the absorbing material has a relatively 
large distance from the GRB. 

We do not detect N V in our spectra but provide an upper 
limit of logA^ < 13.6 6. NV was detected in only four sight- 
lines towards GRBs (Pr ochaska et alJl2008h and likely traces the 
immediate environment of the GRB as it has a high ionization 
potential and requires a strong radiation field. Our upper limit is 
lower than the column densities for those GRB sightlines where 
N V could be detected which again implies that the absorbing 
gas probed by our spectra are most likely not close to the GRB 
itself. Furthermore, we do not detect any fine structure lines and 
only derive an upper limit on Si II* of log N/cnT^ < 12.76. Fine 
structure lines are assumed to be pumpe d by the UV radiation 
field of the GRB (IVreeswiik et al.ll2007h which also indicates 
that the gas is likely very far from the GRB itself. Si II* requires 
a less strong r adiation field and has been detected in Lyman 
break galaxies (iPettini et al.ll2002h where the UV radiation from 
young stars provides the necessary radiation field. 

6. Discussion and conclusions 

GRB 060526 had a relatively bright afterglow that allowed us 
to obtain a solid dataset, both photometrically and spectroscop- 
ically. We achieved a dense light curve coverage over several 
days which allowed a detailed study of the afterglow properties, 
and obtained a series of low resolution but high signal-to-noise 
spectra to study the host environment. 



a elements are produced in massive, metal-poor stars through the a 
process and include elements with integer multiples of the He nucleus 
mass such as O, Si, S, Ca, Mg and Ti. 
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The optical light curve can be fitted with a double smoothly- 
broken power law with a breaks at thi - 0.090 + 0.005 and f/, = 
2.216 + 0.049 days, and decay slopes of apiateau - 0.288+0.026, 
ai = 0.971 ± 0.008, and ai = 2.524 + 0.052. The dense sam- 
pling of especially the 7?c-band light curve also reveals addi- 
tional variability on top of the power laws. These features could 
be explained either by extended activity of the central engine or 
through interactions of the shock with the interstellar medium. 
For the case that the variability arises from external shocks, sev- 
eral mechanisms have been considered. In GRB 021004, both 
density variations of the external medium into which the GRB 
jet plows and angular inhomogenities of the jet surface were con- 
sidered (Nakaretal. 2003). However. iNakar & Grano3 (l2007h 
show that density variations would cause much smaller fluctu- 
ations than those observed in GRB afterglows and can there- 
fore be ruled out. Another possibility is the injection of ad- 
ditional energy into the shock by slower shells that catch up 
with the shocked region as it decelerat es, this model was used 
succe ssfully to describe GRB 021004 (Ide Ugarte Postigo et alJ 
|2005|) and also works better than two other models (double jet 
and density fluctuatio ns) to describe the highly complex light 
curve of GRB 030329 (iHuan^ et al.ll2006l) . Thus, variability can 
give either information on the medium surrounding the GRB or 
on the activity of the central engine. A more intriguing pos- 
sibility is that the flares may be emitted from another region 
closer to the central engine, resulting from late internal shocks. 
Powerful X-ray flares that are attributed to late central engine 
activ ity have been observed in about 50% of al l Swift GRBs 
(e.gjBurrows et all llool IChincarini at all 120071: iKrimm et all 
|2007|) . and strong optical/NIR flaring contemporaneous with the 
GRB pro mpt emission may also occur ( Vestran^ et aL 2005i 
iBlake et al. 2005; Vestrand et al.l2 006; Rac usin et al.l2'008l) . thus 
making optical flares from late central engine ac tivity an inter- 
esting prospect (Kann 2007; Mal esani et al.ll2067h . 

Indeed, iDai et al.l (120071) have suggested that the optical 
variability of the afterglow of GRB 060526 is due to flares 
from late internal shocks (the very early rapid optical variabil- 
ity we present here is very probably due to central engine ac- 
tivity, as it is seen con temporaneously in gamma and X-rays). 
iKhamitov et aP (l2007h . on the other hand, conclude that the 
short timescale of the variabilities requires the jet to be non- 
relativistic already at ~ 1 day and could then be explained by 
external density fluctuations. Our analysis lends tentative sup- 
port to the notion of flares from internal shocks, finding de- 
cay slopes for two flares that exceed what should be possible 
from external shocks. But we caution that the errors of these 
fits are large due to a low amount of data in the decaying parts. 
Furthermore, globally, a model using refreshed external shocks 
is able to account for the light curve variations, although mi- 
crovariability remains. T his creates the i ntriguing possibility of 
reverberation effects (see ' Vestrand et alll2006l for a case of re- 
verberation between gamma-rays and optical emission). Short 
flares in the X-ray or optical bands signal internal shocks from 
long-term central engine activity, and when these shells catch 
up with the forward shock front, they re-energise the external 
forward shock. The detection of such behaviour would probably 
require dense multi-band observations of a bright afterglow to 
search for SED changes at high time resolution combined with 
detailed modelling of the data. This way, one could discern be- 
tween internal shocks (which are expected to have a different 
spectral index from the forward shock afterglow) and refreshed, 
external shocks (which are achromatic). Our data set of the af- 
terglow of GRB 060526 does not allow us such a detailed de- 
composition. 



From the analysis of our low resolution spectra with different 
resolutions, we detect a LLS and a number of metal absorption 
lines that all lie at a redshift of z = 3.221. The low resolution 
only allows us to derive column densities from measuring the 
EWs of the absorption lines and adopt a MISC-CoG analysis 
where we exclude the most saturated as well as blended transi- 
tions. We find a best fit for the Doppler parameter of = 39 ± 3 
km s ' and most of the ions used for the fit lie on the linear part 
of the CoG which allows a relatively reliable determination of 
the column densities. The relative abundances of different met- 
als in the spectra indicate some dust extinction, but an intrinsic 
difference due to enhancement of the production of certain ele- 
ments cannot be excluded. The very low amount of dust detected 
in the afterglow SED may indicate that the latter might be the fa- 
vored possibility. 

The column density of neutral hydrogen is rather low com- 
pared to other GRBs. We derive a metallicity for the host of 
[Fe/H] = -1.09whichis slightly higher than metallicities deter- 
mined from other GRB afterglow spectra. According to the defi- 
nition of QSO absorbers, the host of GRB 060526 is classified as 
a LLS (19 < logA^Hi/cm"^ < 20.3), which seem to have on av- 
erage higher metall icities than damped Lyman a systems (DLA; 
Per oux et al.ll2007 l) and a steeper evolution towards lower red- 
shifts. Around redshift 3, however, the metallicities of both sam- 
ples are within the same range. Also, GRB hosts s how a trend 
towards increasing metallicity with lower redshifts ( Fvnbo et al.l 
2006; Savaglio 2006). Taking into account that most of the sam- 
ple used only low-resolution spectra to derive the metallicity 
(which only gives lower limits for the column densities and the 
metallicity) this evolution might, however, not be as pronounced 
as for DLAs and LLS. This might imply that the enrichment of 
the ISM in the early universe had taken place at earlier times 
than assumed. The absorbing material along the line-of-sight is 
mostly in the neutral state, as usually observed for long GRB- 
DLAs, while sightlines with lower logA^ni often contain more 
ionized material. This might either imply that we have a very 
small host galaxy or that the GRB is placed somewhere in the 
outskirts of its host. In general, the low ionization points to a 
relatively large distance of the absorbing material from the GRB 
itself. 

There is no underlying host galaxy of GRB 060526 detected 
down to a deep limit of 28.5 mag (in F775W AB) in HST/ACS 
data. At that redshift, this means the host has an absolute mag- 
nitude M(1500A) > -18.3 mag, fainter than an 0.5 L* galaxy at 
that redshift. Long GRBs have been found to occur in actively 
star-forming galaxies and star formation is assume d to shift to- 
wards smaller and fainter galaxies over time (e.g. ICowie et al.l 
1996) while massive galaxies prove to be rather unchanged 
through out the histor y of the universe (e.g. lAbraham et ani999l : 
Heavens et alj |2004|) . One would therefore expect that GRB 
hosts should also have h igher luminosities towards higher red- 
shifts. Fvnbo et"an (1200 8) concluded that the observed metal- 
licity distribution of GRB hosts (as well as QSO absorbers) at 
z ^ 3 can be explained by the luminosity function of galaxies 
at that redshift and assuming a luminosity -metalhcity relation 
as derived for other high-redshift samples dLedoux et al.ll2006l ; 
lErb et al.1 '2006). The non-detection of the host of GRB 060526 
down to deep limits, however, would not support this suggested 
evolution. The data neither strongly support nor allow us to rule 
out that the GRB is associated with the nearby galaxy. While the 
offset would be very large compared to typical long GRB offsets, 
it is possible the burst occurred in a locally dense star-forming 
region which is not detected even in our very deep imaging. 
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Table 1. Broad band observations of GRB 060526, times given are the midpoint of the observations. Data have not been corrected for Galactic 
extinction. All data except the Keck/LRIS g' and the HST F115'W data points (in AB magnitudes) are in Vega magnitudes. 
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Notes. UVOT is the 30 cm Uhra Violet and Optical Telescope onboard the Swift satelUte. 
Maidanak is the 1.5m telescope of the Maidanak observatory in Uzbekistan. 

RTT150 is the 1.5m Russian-Turkish Telescope at TUB iTAK National Observatory on Mount Bakyrlytepe, Antalya, Turkey. 
MOSCA is the MOSaic CAmera on the 2.5m Nordic Optical Telescope, La Palma, Canary Islands, Spain. 

Keck/LRIS is the Low Resolution Imaging Spectrograph on the 10m Keck I telescope, Mauna Kea, Hawaii, United States of America. 
DOLORES is the Device Optimized for the LOw RESolution detector on the 3.6m TNG (Telescopio Nazionale Galileo) telescope on La Palma, 
Canary Islands, Spain. 

MIRO is the 1. 2m telescope of the Mt. Abu Infrared Observatory, India. 

BFOSC is the Bologna Faint Object Spectrograph & Camera at the G. D. Cassini 152 cm telescope of the Bologna University, Loiano, Italy. 

DFOSC is the Danish Faint Object Spectrograph and Camera on the Danish 1.54m telescope on La Silla, Chile. 

Watcher is the 0.3m Watcher robotic telescope at Boyden Observatory, South Africa. 

Shajn is the 2.6m Shajn telescope of the Crimean Astrophysical Observatory (CrAO), Ukraine. 

TLS 1.34m is the 1.34m Schmidt telescope of the Thiiringer Landessternwarte Tautenburg, Germany. 

ALFOSC is the Andalucia Faint Object Spectrograph and Camera on the 2.5m Nordic Optical Telescope, La Palma, Canary Islands, Spain. 

FORS2 is the FOcal Reducer and low dispersion Spectrograph 2 on the 8.2m Very Large Telescope, Paranal Observatory, Chile. 

ROTSE is the 0.3m Robotic Optical Transient Search Experiment III-B telescope at the H.E.S.S. site, Mt. Gamsberg, Namibia. 

ANDICAM is the A Novel Double-Imaging CAMera detector on the 1.3m Small and Moderate Aperture Research Telescope System (SMARTS) 

telescope at the Cerro Tololo Interamerican Observatory in Chile. 

PAIRITEL is the 1.3m Peters Automatic InfraRed Imaging TELescope on Mt. Hopkins, Arizona, United States of America. 
HST ACS is the Hubble Space Telescope equipped with the Advanced Camera for Surveys. 



